Effective t-J Hamiltonian for the Copper Oxides. 
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Abstract 

Starting from the Emery model, which is assumed to describe the cop- 
per oxygen planes, and including direct oxygen hopping matrix elements, t pp , 
we have been able to derive the effective t — J Hamiltonian for the copper 
orbitals using the Linked Cluster Expansion Method up to fourth order in 
tpd, the hybridization matrix element. The spin depend part of the effective 
Hamiltonian is composed of two contributions: the superexchange interaction 
and another one of RKKY-type. The effective parameters (t and J) depend 
on doping, 8. This Hamiltonian can be used to study the magnetic properties 
of the High-Tc materials vs. 8. 
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The study of High-Tc materials continue to attract attention both from the the- 
oretical as from the experimental point of view/ 1 ' There is general agreement that 
the basic structure of these materials is composed of the CuOi planes. There is 
evidence from both band - structure calculations and spectroscopy investigations 
that the generic properties of the CuOi planes can be described only in terms of the 
Cu3d x 2_ y 2 and the Op x [ y ) orbitals. These considerations lead us to the three - band 
Hubbard or Emery model. This is given by/ 2 ' 

H = (e d -fi)^2 d\ a d- la + (e p - /i) E pLpi<? + u dJ2 4t d it4i d h 

\a \a i 

+t PP E (-1) % pU„ + ^ E (-l) M K(pid ia + H.c.) , (1) 

<1,1'><7 <l,i>a 

where t p d is the hopping matrix element between the Cu and O orbitals, Ud is the 
Coulomb repulsion of two holes in a Cu d x i_ y i orbital, A = e p — e^, or charge trans- 
fer gap (CT), between the p and d orbitals, in the hole picture and t pp is the direct 
oxygen hopping matrix element. These are the most important parameters of the 
model. These parameters are usually taken from LDA calculations® to be : A % 
3.5 eV, tpd « 1.35 eV, t pp « 0.65 eV, and Ud ~ 9eV. In Eq. (1), the indices 1 and i 
refer to the oxygen and copper sites, respectively. \i is the chemical potential. The 
phase factors and M\\ are either or l/ 4 ' 

The purpose of this letter is to give a derivation of the t — J® starting from the 
Emery model of Eq. (1). The t — J has excitations which correspond to the low - 
energy excitations of the Emery model/ 1 ' 

To start our derivation, we Fourier analyze the oxygen orbitals. This gives: 



H = (e d - n) E 4<? d i<? + E e kvcl^Ckva + U d Yl 4t d it4i d H + H ' > ( 2 ) 

ia \lvg i 

where H' is given by 

^ £ M e ~*'M,<W + (3) 
v N ikwT 

and 

e kv = (e p -n) +4t pp (-iysin^-sin^- ; 7kl , = ^(sin^ + {-If sin^) 1 ' 2 , (4) 
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In Eq. (3), N is the number of copper sites. This model is the three band Ander- 
son lattice model with the explicit oxygen band, e^J 4 ' In contrast to the heavy - 
fermion systems, the oxygen bandwidth, 4t pp , has about the magnitude of t pd and 
the oxygen bands are lightly filled for 8 — > 0. 

It can be shown that it is a good approximation to reduce the two oxygen bands 
to a single oxygen band/ 6 ' This is done by means of a rotation given by Ck CT = 
~ k ^k"k.T- where has the following form 



M k = 



7ki 

7k2 



7k2 

-7ki 



where 



7k = 7ki + 7k2 = 4 [sin {—) + sin { — 



(5) 



Then, we end up with the following two band Hubbard Hamiltonian 



H = (e d - fi)Y, 4a d ™ + e k^L^k CT + U d 4t d it d li d il + H ' j ( 6 ) 

\a k<7 i 

where H' is given by 

^'^T=E7k(^ k - R i4^ + ^ C .) (7) 

and 

e k = (e p - /i) + At pp \sin-^-sin^-\ 7je = 2{sin 2 -^ + sin 2 ^-) 1/2 , (8) 

where |F| means absolute value of F. Let us mention that we are going to work 
with one oxygen band for simplicity. The case of two oxygen bands can be treated 
similarly, but the final expressions become more involved. We are trying to keep the 
algebra as simple as possible. 

Next, we want to obtain an effective Hamiltonian for the coupling between the 
copper atoms. For this purpose we use the Linked Cluster Expansion Method/ 7 - 1 
This theorem tells us that the thermodynamic potential, J7, can be calculated from 
the following expression: 

-1 oo 

n = a, - -J2 u n , (9) 
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where fl is the unperturbed thermodynamic potential, (3 is the inverse of the tem- 
perature and U n is given by 



_ (-iy 



n 



dr x ... \ dr n <T T H'(r 1 )...H'(r n ) > c , 



(10) 



and T T means time ordering. The averages are performed with the oxygen orbitals 
and c means different conected diagrams. 



Since, 



<^: +1) (r)>^<^ n+1 \r)>=0 , (11) 

only even terms of Eq. (10) contribute. Thus, in second order in t p d, we have, with 
H 2 = U 2 , 



Ho = 



E [ d ^[ dT ^' [^ (k ' Ri+k '' R ' i ' ) ^^L(^)4 CT (^)>^(^)^ CT '(^) 



2/3N TT Jo Jo 



i'kV 



+ ^ (k - Rrk '- R 'i' ) (T T 4 CT (r 1 )0 kCT (r 2 ))d iCT (r 1 )4 CT ,(r 2 ) 
+ e i(k ^- k '^'\T T ^(r 1 )^(75 ! )>4 T (r 1 )d iV (7 5 ) 
+ e 4(k ' Ri+k '' R i' ) (T^ kCT (r 1 )^ kCT (r 2 ))4(r 1 )4 CT ,(r 2 



But, 



{<ff(fj) c ={<M>)c = o , 

{TAM^'M)) = ^S^iT^ljr^ir,)) ,etc 
Thus, Eq. (11) reduces to: 



(12) 

(13) 
(14) 



Ho 



i'kV 



+e ik.R ;; , < T^inmjr,) > c d\ a {n)d- w {ri 



Since, 



G°(k, t) = - < T T kCT (r 1 )4 <7 (r 2 ) > c • 
So, Eq. (11) is transformed to, 

h 2 = Y: t eff (n)(did- l+n „ + H.c) , 



(15) 



(16) 



(17) 
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where, the effective hopping matrix element, t e ff(H) is given by 



W/(R)^E7£cos(£.i?) 



(n F (e d + U d ) - n F (e^)) 



(n F (e d )-n F (e n )) 



(18) 



(^ - £d ~ U d ) 

with n F (e) being the Fermi - Dirac distribution function and (R) means nearest 
neighbor atoms (n.n.). 



Next, let us evaluate the fourth order term, H4 = U4. It is given by 
H 4 = -— •••/ <T T H (n)---H (r 4 )> c • 



4/? Jo Jo 

After a straighforward calculation, we arrive to the following result 
ti 



(19) 



4iV 2 



E (7k7k') 2 ^(k,k';/3) 



ikp, 
i'k'p' ' 



P<7 1 CT p' CT 



where 



e - 4 (kR ;; ,+k R p , p ) | | 

1(7 p<7 P CT 1 CT 



i(k,k';/3) = ^ G c(k,iA: n )G°(k',iA; n )^(^V 



and 



G°(k,iA: n ) = 



(ik n - ek) 



+ 



2(i&„ - e d ) 2(iA; ri - q - C/ d ) 
where we have taken the atomic limit for the copper sites. ^ 

By performing the Matsubara summation, we obtain 

A(k,k';P) = A 1 (k,k';P)+A 2 (k,k';!3) , 



(20) 

(21) 

(22) 
(23) 



(24) 



with 



ii(k,k';/?) 



2(e k -Q)(e k ' -e d ) L(e k ' ~ f 



n F (e rf + U d 



2(e d + U d -6 k )(e d + U d -6 k , 



(ek - Q - U d ) 



+ 



(25) 



i 2 (k,k';/?) 



2(fk - 6 k ' 



e k - €d) 



+ 



+ 



(e k - e d - U d f ' (6k - e d ) (e k - e d - U d ) J ' (26) 

In order to obtain a cos(k — k ) • R dependence from Eq. (20), we need to make 
the following identifications 

i = p ; p = i (first term), 
i = p ; p = i (second term), 
i = p ; p = i (third term, etc.) 
Then, we arrive to the following result 
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H 4 = o E J ( R M-4 CT ' d i+R CT ,f 4+R CT 



i<R> 



(27) 



where 



3t 

J ( R ) = ^E i ( k ' k '^)(7k7 k 2 cos(k-k')-R 
kk' 

On the other hand, the spin summation can be rewritten as 



(28) 



E^L'V^ = 2 n;nj + 2Si ' Sj ' (29) 

(7(7 

So, H4 can be expressed in a compact form as: 



H 4 = ^ J(R) Sl • S ;+R + -n ; n ;+R • (30) 

i<R> V 4 J 

From Eq. (30), we can conclude that the fourth order contribution to the effec- 
tive Hamiltonian has the form of spin dependent part J(R), as in the t — J model. 
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This J(R) term has two contributions (See Eqs. (24), (25) and (26)). If we make 
the following approximation 



eg - e d « A , (31) 
we see that A m (k, k';(3), m = 1, 2 is given by 

ii(k,k';/?) « D x , (32) 

i 2 (k,k';/?) «£ 2? , (33) 

where D\ and D2 are constants, with D\ > 0. 

J(R) has two contributions: 1) a superexchange type of interaction due to 
ii(k, k'; (3); {9) and 2) a RKKY-type of interaction ( 6 ) (i 2 (k, k'; (3)) which is strongly 
dependent upon doping rate, 5 and has an oscillatory behavior at very long distance 
between the copper atoms. It is this doping dependency of this J(R) which makes 
it interesting. It is possible to study the magnetic properties as function of doping, 
as it was done in ref. [6]. In particular, the Jrkky interaction might lead to a 
frustation of the long - range AF ordering with doping (small 5). In addition to the 
J(R) term, we find a local repulsion given by the Ud term. 

At this point, we have the effective Hamiltonian, H e ff, on the copper sites, given 

by 

H eff = (e d -fi) J2 d lA° + H 2 + H 4 + UdY,4t d iAi d H • (34) 

ia i 

In the large Ud limit, we can enforce no double occupancy at half filling, by using 
perturbation theory on virtual states. ( 10 ) The result is 

Heff = (Q-^E^L^ + ^E d\+Radia 
\a i(R)<T 

+ ( Ji + J 2) E S; • S; + j^ + ( Ji — J 2) E n i+R n i + (three site terms) , (35) 
i<R> i<R> 

where 

t = t eff (R) ; Ji = J(R) ; J 2 = — , (36) 

dia = d- la (l - n;_ CT ) • (37) 

In conclusion, we have been able to derive the t — J Hamiltonian starting from 
the three band Hubbard or Emery model. Our effective parameter Ji contains the 
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superexchange interaction of the t — J derived by Zhang and Rice/ 5 ) in the limite of 
t pp — y 0. However, t pp is a key element in the description of the physics of High-Tc 
superconductors and it should not be taken equal to zero. The effective parameters 
depend on doping by means of the Fermi level. The spin dependent part of the ef- 
fective Hamiltonian is composed of two contributions, one mainly antiferromagnetic 
given by A\, and a second one given by A2, which is of the RKKY - type. To our 
knowledge, it is the first time that the latest contribution has been derived simulta- 
neously with the superexchange contribution. There is a calculation by Eskes and 
Jefferson^ 11 ) where they have included t pp and then they use fifth order perturbation 
theory. However, they do not obtain the RKKY - type interaction. The presence 
of these two contributions are interesting for studying the evolution of the magnetic 
properties upon doping. 
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